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A new investigation is done of the possibility of binding the "dodecatoplet", a system of six top quarks and six 
top antiquaries, (i 6 ^ 6 ), using the Yukawa potential mediated by Higgs exchange. A simple variational method 
gives a upper bound close to that recently estimated in a mean-field calculation. It is supplemented by a lower 
bound provided by identities among the Hamiltonians describing the system and its subsystems. 
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O ' I. INTRODUCTION 

o : 

The Higgs boson, responsible for generating the masses of fermions in the standard model, couples more strongly to the heavy 
O . quarks than to the light ones. The question has been raised whether the attraction mediated by Higgs exchange could produce 
new type of bound states [1—5]. Higgs exchange corresponds to the Yukawa type of potential 

m i -a H v(r) = _ Qg cx P(-H . (1) 

(N . r 

, 1 As discussed in the literature [1-3], the coupling of Higgs to the top quark should be about g t ~ 1, with au = gl /(47r). We 

■ thus consider an — 1/(4tt) as a benchmark value for estimating the spectrum of the potential v, but we shall also study how 
O-f the eigenenergies behave as a function of an- As for /i, the Higgs mass, it expected to be is of the order of a hundred or a few 
jV- hundreds of GeV. 

If a bound state of several quarks and antiquarks occurs due primarily to the above potential, the precise determination of its 
i-C . mass in a regime of strong binding should incorporate relativistic effects, strong forces, M^-exchange inducing it <-> bb mixing, 
etc. [3, 5]. However, in a regime of weak binding the system remains non-relativistic. Thus for a preliminary investigation of 
the existence of new type of bound states due to Higgs exchanges, it is sufficient to rely on the Hamiltonian [1, 2, 4] 

N 2 

Hn = J2 ^ ~ aH w ( r «) ' (2) 
r***'" ■ » i<j 

and examine its spectral properties. If N, the number of constituents, does not exceed 6 top quarks and 6 top antiquarks, 
the colour and spin degrees of freedom can endorse the constraints of antisymmetrisation, and for the orbital variables, the 
Hamiltonian (2) can be considered as acting on effective bosons. This is why the attention has been focused on the (i 6 F) 
<^> system, which can be named "dodecatoplet", by analogy which the late pentaquark. 
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II. UPPER BOUNDS 
A. Scaling 



It is well-known that under the (xri — > r t transformation for all positions, and thus for the interquark distances r^, the level 
energies of Hn, and in particular its ground-state, scale as 

tp ( \ M (r <\ n m t a H 

E N {rn t ,a H ,Li) = — e N (&) , &= . (3) 

m t fx 

Thus one is dealing with a one-parameter problem. For ajj = g , |/(47r), m t = 172.6 GeV, a reasonable Higgs mass near 
fi = 140 GeV, one gets G ~ 0.1 as the order of magnitude of the dimensionless coupling. 



B. The two-body case 



For two bosons, in a Yukawa potential, after scaling the variables and removing the centre of mass, the problem is reduced to 

h 2 = -A-Gexp(-r)/r , (4) 
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FIG. 1: ejv/(A r — 1), where ejv is the A^-body energy problem as a function of GN/2, where G is the coupling. Solid curve: exact 2-body 
energy e2 which is also a lower bound to ejv/ (-/V — 1) for N > 2. Dashed curve: Gaussian variational approximation £n/(N — 1), valid for 
any N. o:variational hyperscalar approximation £3/2 for the 3-body case. •: hyperscalar approximation £12/11 for the 12-body case. 



which starts supporting a bound state for G > G2 — 1.68, as shown in the classic paper by Blatt and Jackson [6]. With G ~ 0.1 
we are thus far from a Yukawa binding of (tt). The ground-state energy €2 is drawn in Fig. 1 as a function of G. Also shown is 
the variational upper bound e 2 = xmn a [t(a) — Gp(a)] corresponding to a single normalised Gaussian, tl> a ( r ) ^ exp(— ar 2 /2), 
whose range parameter a is optimised. The relevant expectation values, involving the complementary error function, are 



t(a) 



■A) -I 



p{a) = (v(r)) 



2y/a 



exp 



1 

17i 



erfc 



2^ 



(5) 



This trial function demonstrates binding for G > 2.71 only, to be compared to the exact G > G2 — 1.68. This gives an idea 
of the validity of a simple variational approximation in a regime of weak binding. A function with better behaviour at large r 
would of course do much better. 



C. Simple variational upper bound for the iV-body system 

If one describes the relative motion using (N — 1) Jacobi variables 

2r 3 - n - r 2 

xx=r 2 -r 1 , x 2 = -j= (6) 

from the rescaled positions, and removes the centre-of-mass motion, the A^-body Hamiltonian reduces to 

w-i 

Y^X G^rir.J. (7) 

i=i i<j 

If the above Gaussian is generalised as Ilta 1 ^a{ x i) anc l taken as trial function, the variational energy reads 

e N = min [(N - l)t(a) - A^(A^ - l)Gp(a)/2] , (8) 

a 

which is easily estimated, and becomes a better and better approximation as increases. Note in (8) the obvious relation 

ejv(G) = (N - 1) i 2 (NG/2) , (9) 

which, as seen below, becomes a constraining inequality if the approximate energies are replaced by the exact ones. It also 
indicates that if the number N of bosons is large enough, binding is achieved however small is the coupling G. This possibility 
of binding systems whose subsystems are unbound, first pointed out by Thomas [7], is nowadays refereed to as "Borromean 
binding" . For references, see, e.g., [8]. 

Note also that the variational energies obey the same scaling laws, virial theorem, etc., as the exact ones. This is stressed 
in [4], and can be traced back up to very early papers dealing with quantum systems [9, 10]. 
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D. Hyperscalar upper bound for iV-body systems 
The above trial wave function J\ i ^ a (xi) is a particular function, namely Gaussian, of the hyperradius r > defined as 



■ x 



N-l ■ 



(10) 



It occurs when the relative distances Xi are considered as a unique vector in a 3 (N — 1) dimensional space and polar coordinates 
{r, ft} are introduced there to solve the Schrodinger equation. See, e.g., [11]. As the potential ^ vfaj ) is usually non central 
in this space, the Schrodinger equation in these coordinates is expressed as an infinite set of coupled radial equations. However, 
for a system of bosons, the ground-state is well described in the (variational!) approximation consisting of retaining the lowest 
partial wave: this is the hyperscalar approximation, which reads 



u"(r) + 



L(L + 1) N(N-l) 
e- „ + G — — J -v 00 (r) 



u{r) = 



(11) 



with e < e since the Gaussian approximation is a particular ansatz for u(r). Here, L — 3(N — 2)/2 is an effective orbital 
momentum, and the hypercentral projection reads 



voo{r) 



,tt/2 



cos 2 9 sin" 9v(r cos 9) d9 



f w/2 

/ cos 2 9 sin"0d6> 
Jo 



(12) 



where n = 3N — 7. There is no difficulty to evaluate woo analytically 1 , but for large N, the analytic expression contains 
cancelling large positive and large negative terms, and a direct numerical estimate of (12) might lead to better accuracy and 
stability in the computation. 



E. Results 



For (unphysical) G = 4 and N = 12, one gets e ~ -1096 vs. e ~ -1117, which translates into (H)/[N(N - l)a 2 H m) ~ 
—0.0480 for n/[(N — l)majj] = 0.022 in the notation of the authors of Ref. [4], very much compatible with their estimate 
based of a variational function JTcxp(— ari) and the result of a Hartree-Fock calculation. The hyperscalar approximation is 
slightly below, i.e., better: the quality of the trial wave-functions in [4] is partly compensated by the lack of removal of the 
centre-of-mass motion. 

The variational energies 1(G) and e(G) in the physical domain (small G) are shown in Fig. ??. A system of 12 bosons starts 
supporting a bound state for G > 0.45 from the Gaussian approximation and the hyperspherical equation. The hyperscalar 
approximation is very close to the single-Gaussian approximation for N = 12, while for the N = 3 case, also shown in Fig. ??, 
it is appreciably better. The estimate of [4], if translated in our notation (this is the hh < 0.19 of their figure, corresponds to 
G > 0.478. Note that Pacheco et al. [12, 13] considered previously "self-Yukawian" boson systems, using a self-consistent 
Hartree method similar to that of [4] and obtained binding for a slightly better G > 0.438. For qlh = l/(47r), this corresponds 
to a maximal Higgs mass of about 31 GeV. 



III. LOWER BOUND FOR THE GROUND-STATE ENERGY 



Even though the above hyperspherical approximation and the central-field method of Refs. [4, 12] are known to be very 
good approximations to the exact energy, the doubt could remain that the existence of bound states has been missed, or that the 
absolute value of the binding energy has been underestimated. It is then desirable to derive a lower bound on the ground-state 
energy. 

A basic tool consists of splitting the Hamiltonian into pieces, say 

H = A + B + ■ ■ ■ => E(H) > E(A) + E(B) + ■■■ , (13) 



1 / voo( r ) r2L+2 exp(— ar 2 ) dr/ J r 2L + 2 exp(— ar 2 ) dr = p(a) provides a good cross-check. 
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in an obvious notation where E(H) is the ground-state energy of H . Saturation is obtained if A, B, etc., reach their minimum 
simultaneously. 

The simplest application to our problem is done by rewriting the TV-body Hamiltonian of (2) as a sum of 2-body ones 

H N [m t ,a H ] = j^-j ^2 H 2 j [ m t,(N-l)a H ] , (14) 

i<j 



and leads to 



or in the rescaled variables 



N 

E N [m u a H ] > — E 2 [m tl {N-l)a H ] , (15) 



N 

e N (G) > -e 2 [(N - 1)G] , (16) 

but this is not very accurate because the energy of each two-body subsystem is bounded by its rest energy, although each pair 
has an overall motion within the whole TV-body system. The remedy [13-15] consists of writing identities among the intrinsic 
Hamiltonians 

ft _ n (Pi + ---+Pn) 2 nr . 



namely 



This leads to the improved inequality 



H N (m t ,a H ) = ^J2 H * j ( mt ' n °h/ 2 ) • (18) 



N 



E N {m, a H ) > {N - 1) E 2 (m, Na H /2) , (19) 
which is always better [15] than the previous inequality (15). In the rescaled variables, it reads 

e N (G) > e N (G) = (N - l)e 2 (/VG/2) . (20) 

For the above numerical example with G = 4, we obtain e 12 ~ —1335, to be compared with the upper bound ei 2 ~ —1117. 
As seen in Fig. 1, the window is rather narrow between the variational upper bound and the Hall-Post lower bound. 

But we are here in a regime of deep binding that would require drastic relativistic corrections. A more important consequence 
of (18) is that Hn hardly supports bound states if each H 2 is positive. This means that binding requires G > Gn with 

G N > |p G 2 ~ 0.28 , (21) 
which corresponds (again for an = 1/ (47r) to a minimal Higgs mass which is certainly below 49 Ge V. 



IV. SUMMARY 



We have revisited the problem of binding 6 top quarks and 6 top antiquarks with a Yukawa potential mediated by Higgs 
exchange. The problem depends only of the dimensionless coupling G = m t aH / ' H- By combining variational estimates and 
Hall-Post type of lower bound, the 1 2-body problem is under control. In particular, it is found that the minimal strength to bind 
the system belongs to the interval 

0.28 < G^ 2 < 0.43 . (22) 

With a Higgs coupling taken as g t = 1, corresponding to a H = l/(47r), the maximal Higgs mass leading to this binding should 
be 



31GeV < M max < 49GeV . 



(23) 
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With the current ideas on the Higgs mass and its coupling, the existence of such dodecatoplet seems very unlikely. Perhaps 
higher systems can be bound. For 3 He atoms, the interaction potential contains attractive parts, but is unable to bind the dimer, 
unlike the 4 Hc case, where the dimer is bound, due to the heavier mass of the constituents. It was found, however, that binding 
becomes possible for N > 35 atoms of 3 He [16, 17]. Of course, to study (t n i m ) systems with n or m larger than 6, one should 
account for the antisymmetrisation effects in building the orbital wave function. 
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